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Zone methods and the fermion sign problem 
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We review a recently proposed approach to the problem of alternating signs for fermionic many body Monte 
Carlo simulations in finite temperature simulations. We derive an estimate for fermion wandering lengths and 
introduce the notion of permutation zones, special regions of the lattice where identical fermions may interchange 
and outside of which they may not. Using successively larger permutation zones, one can extrapolate to obtain 
thermodynamic observables in regimes where direct simulation is impossible. 



1. Introduction 

The zone method approach to the fermion sign 
problem is based on the observation that in many 
finite temperature simulations fermion permuta- 
tions are short ranged This can be true 
even for systems with massless modes and long 
distance correlations, as we demonstrate with 
an explicit example. In this proceedings article 
we derive an estimate for the finite temperature 
fermion wandering length and discuss the main 
features of the zone method. 

2. Worldlines 

We begin with a brief review of the worldline 
formalism We introduce ideas in one spa- 

tial dimension before moving on to higher dimen- 
sions. Let us consider a system with one species 
of fermion on a periodic chain with L sites, where 
L is even. Aside from an additive constant, the 
general Hamiltonian can be written as 

H = [aj^ifl,, + a|ai+ij +^ Ciojai. (1) 
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Following d we break the Hamiltonian into two 
parts. He and Ho, 
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*The author thanks Shailesh Chandrasekharan, Hans Ev- 
ertz, Pieter Maris, and Nathan Salwen for discussions and 
comments. This work benefitted from the resources of the 
North Carolina Supercomputer Center. 



We note that H 
write 



He + Ho- For large N, we can 
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where 



Sc/o = exp(-;^i7o/o)- 



(2) 



(3) 



Inserting a complete set of states at each step, we 
can write Tr [exp{—f3H)] as 

^ {zn\So\z2N-i) ■■■{zi\Sc\zo) ■ (4) 

The worldline trajectory of each of the fermions 
can now be traced from Euclidean time i = to 
time t = p. Since we are computing a thermal 
trace, the worldlines from t = to t = f3 define a 
permutation of identical fermions. Even permu- 
tations carry a fermion sign of -1-1 while odd per- 
mutations carry sign —1. The generalization to 
higher dimensions is straightforward. In two di- 
mensions, for example, exp(— /3i/) takes the form 



exp 
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(5) 



3. Wandering length 



In one spatial dimension, the Pauli exclusion 
principle inhibits fermion permutations except in 
cases where the fermions wrap around the lattice 
boundary. For the remainder of our discussion, 
therefore, we consider systems with two or more 
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dimensions. The first question we address is how 
far fermion worldhnes can wander from start time 
t = to end time t = (3. We can put an upper 
bound on this wandering distance by considering 
the special case with no on-site potential and only 
nearest neighbor hopping. 

Let us consider motion in the x-direction. For 
each factor of S^S^ in a given fermion may re- 
main at the same x value, move one lattice space 
to the left, or move one lattice space to the right. 
If h is the hopping parameter, then for large N 
the relative weights for these possibilities are ap- 
proximately 1 for remaining at the same x value, 
(3hN^^ for one move to the left, and j3hN~^ for 
one move to the right. In (||) we see that there 
are A'' factors of S'^S^. Therefore for a typical 
worldline configuration at low filling fraction, /, 
we expect ~ (3h hops to the left and ~ ph hops to 
the right. For non- negligible / some of the hops 
are forbidden by the exclusion principle. Assum- 
ing random filling we expect [1 — f)(3h hops to 
the left and ~ (1 — /)/3/i hops to the right. 

The net displacement is equivalent to a ran- 
dom walk with 2(1 — f)f3h steps. The expected 
wandering length, A, is therefore given by 

A = ^2{1 -f)(3h. (6) 

In cases with on-site potentials, fermion hop- 
ping is dampened by differences in potential en- 
ergy. Hence the estimate serves as an upper 
bound for the general case. We have checked the 
upper bound numerically using simulation data 
generated by many different lattice Hamiltonians 
with and without on-site potentials. 

4. Permutation Zone Method 

Let W be the logarithm of the partition func- 
tion, 

W = log {Tr [exp(-/3i?)]} . (7) 

Let us partition the spatial lattice, F, into zones 
Zi, Z2, Zk such that the spatial dimensions of 
each zone are much greater than A. For nota- 
tional convenience we define Zq = For any 
-R C F, let Wr be the logarithm of a restricted 
partition function that includes only worldline 



configurations such that any worldline starting at 
t = outside of R returns to the same point at 
t = p. In other words there are no permutations 
for worldlines starting outside of R. We note 
that Wr = W, and is the logarithm of the re- 
stricted partition function with no worldline per- 
mutations at all. Since the zones are much larger 
than the length scale A, the worldline permuta- 
tions in one zone has little or no effect on the 
worldline permutations in another zones. There- 
fore 

Wz,u...z, - Wzou...z,^^ ~ Wz, ~ W^. (8) 
Using a telescoping series, we obtain 

Wr^W^= (W^Zou...z,-M/z„u...z,_J, (9) 

]=i,...,k 

Wv-W^^ {Wz,~W^). (10) 

j=i,...,k 

For translationally invariant systems tiled with 
congruent zones we find 

M/ = VFr«VF0 + ^(VFz, -H^g), (11) 

where |r| / \Zi\ = k, the number of zones. For 
general zone shapes one can imagine partitioning 
the zones themselves into smaller congruent tiles. 
Therefore the result (^0|) should hold for large ar- 
bitrarily shaped zones. For this case we take |r| 
to be the number of nearest neighbor bonds in the 
entire lattice and \Zi \ to be the number of nearest 
neighbor bonds in the zone. We will refer to \Zi\ 
as the zone size of Zi. This is just one choice for 
zone extrapolation. A more precise and compli- 
cated scheme could be devised which takes into 
account the circumscribed volume, number of in- 
cluded lattice points, and other geometric quan- 
tities. 

As an example of the permutation zone 
method, we compute the average energy {E) 
for a free fermion Hamiltonian with only hopping 
on an 8 X 8 lattice. We set the number of time 
steps = 8 and consider values f3h — 1.0, 1.5, 
and 2.0. The corresponding values for A are 
1.0, 1.2, and 1.4 respectively. The Monte Carlo 
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Zone Extrapolations for < E>h (8x8x8) 

PA= 1.0 




Figure 1. Average energy at (ih — 1.0, 1.5, 2.0 for 
the free fermion Hamiltonian on an 8 x 8 lattice. 



updates were performed using our own version 
of the single-cluster loop algorithm ||^. In Fig- 
ure 1 we show data for rectangular zones with 
side dimensions x 0, 1 x 1, 2 x 1, 2 x 2, 3 x 2, 
3x3, 6x6. We also show a least-squares fit 
(not including the smallest zones 0x0 and 1x1) 
assuming linear dependence on zone size as pre- 
dicted in ([l^). We find agreement at the 1% level 
or better when compared with the exact answers 
shown on the far right, which were computed us- 
ing momentum-space decomposition. 

While the physics of the free hopping Hamil- 
tonian is trivial, the computational problems are 
in fact maximally difficult. The severity of the 
sign problem can be measured in terms of the av- 
erage sign, <Sign>, for contributions to the par- 
tition function. For (3h = 1.0, <Sign>-- 0.005; 
for I3h = 1.5, <Sign> - lO'^; and for /3/i = 
2.0, <Sign>'^ 10"^. Direct calculation using 
position-space Monte Carlo is impossible by sev- 
eral orders of magnitude for fih > 1.5. 



is short-ranged and has a maximum range of 
^2(1 — lattice sites, where /3 is the inverse 
temperature, h is the hopping parameter, and / 
is the filling fraction. We have introduced the no- 
tion of permutation zones, special regions of the 
lattice where identical fermions may interchange 
and outside of which they may not. Using succes- 
sively larger permutation zones, one can extrap- 
olate to obtain thermodynamic observables. Ap- 
plications of the zone method to chiral symmetry 
breaking in (2+1) dimensional QED is discussed 
in another article in these proceedings Q. 
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5. Summary 

We have reviewed the zone method approach 
to the fermion sign problem. We have demon- 
strated that the exchange of identical fermions 



